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Modulation-assisted tunneling in laser-fabricated photonic Wannier-Stark ladders 
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We observe Wannier-Stark localization in curved photonic lattices, realized using arrays of evanes¬ 
cent ly coupled optical waveguides. By correctly tuning the strength of inter-site coupling in the 
lattice, we observe that Wannier-Stark states become increasingly localized, and eventually fully 
localized to one site, as the curvature of the lattice is increased. We then demonstrate that tun¬ 
neling can be successfully restored in the lattice by applying a sinusoidal modulation to the lattice 
position, an effect that is a direct analogue of photon-assisted tunneling. This precise tuning of the 
tunneling matrix elements, through laser-fabricated on-site modulations, opens a novel route for the 
creation of gauge fields in photonic lattices. 
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Introduction Quantum matter or light propagating in 
engineered lattices offer versatile platforms for the quan¬ 
tum simulation of new states of matter, such as topolog¬ 
ical phases HHS]. This approach relies on novel technolo¬ 
gies allowing for the tuning of microscopic parameters 
that characterize lattice models of interest, such as on¬ 
site interactions and tunneling matrix elements. In fact, 
generating complex tunneling matrix elements can poten¬ 
tially induce non-trivial gauge structures in the lattice, 
e.g. artificial magnetic fields and non-Abelian gauge po¬ 
tentials offering an interesting route for quantum 

simulation [ 7 ]. In this context, photon-assisted tunneling, 
a powerful method by which tunneling can be controlled 
in lattice systems, has been recently exploited in cold 
gases [SHU] and ion traps da US); this led to the experi¬ 
mental realization of the Hofstadter model miiisj and to 
the detection of the topological Chern number [19] with 
cold atomic gases. 

The photon-assisted-tunneling method relies on two 
main ingredients uniiisHia: (a) an artificial electric field 
generating a large energy offset A between neighboring 
sites, hence inhibiting the bare hopping, and (b) a time- 
modulation of the on-site energy, whose frequency is res¬ 
onant with respect to the static offset uj = A/h; this 
restores the tunneling in an efficient and tunable man¬ 
ner. In this Letter, we experimentally demonstrate the 
realization of modulation-assisted tunneling, an analogue 
of photon-assisted-tunneling, in arrays of coupled optical 
waveguides. 

The transport of light in a system of coupled optical 
waveguides, a photonic lattice, can be described by a 
Schrodinger-like equation. As a result, photonic lattices 
can be used to observe phenomena known from solid state 
physics. In recent years photonic lattices have been used 
to study fundamental solid state phenomena including 
Bloch oscillations ISQlEi], dynamic localisation [22ll24] . 
Bloch-Zener dynamics [25|, and Landau-Zener dynamics 


[26] . These phenomena are each related to the manner in 
which a charged particle behaves in a periodic potential 
and external electric field. In such a system, a static 
electric field destroys the translational symmetry of the 
lattice, and the delocalised Bloch states become localised 
in space. These states are known as Wannier-Stark (W- 
S) states, originally predicted by G. H. Wannier in 1960 

m- 

In the absence of an external electric field, the eigen¬ 
states of an electron in a periodic potential are the Bloch 
states. A static electric field {Sdc) destroys the degen¬ 
eracy of these spatially delocalised Bloch states. In this 
situation, the eigenstates and eigenfunctions are [28l [29| 


^ ; Jn—m 1 ^); 
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where Ju{') is the Bessel function of order n, a is the 
lattice spacing, An is the bandwidth, e is the electronic 
charge and |n) are the Wannier states. The span of the 
first Brillouin zone is —Trja < k < ir/a. In the limit 
2n/e£dcCi 0 there is only one term in Eq. 0 , i.e. 
the eigenstates exactly correspond to the localized Wan¬ 
nier states \m). In fact, in a strong external electric field 
and weak inter-site interaction i.e. 2n/e£dcCi 1, the 
spatial width of W-S state is less than the inter-site sep¬ 
aration, a. In this limit, the W-S states will be localised 
to a single lattice site, indicating that the energy off¬ 
set A = eSdcCi ^ n. generated by the static electric field 
inhibits the bare hopping between neighboring sites. Im¬ 
portantly, when driving the system, the strongly local¬ 
ized electronic states on individual lattice sites can inter¬ 
act with photons and tunnel to the nearest lattice sites. 
This type of tunneling with discrete energy exchange is 
known as photon-assisted tunneling; it has been observed 
in superconducting diodes [30], semiconductor superlat¬ 
tices [31], quantum dots [32], and also with cold atoms 
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FIG. 1: Schematic diagrams of (a) an array of straight op¬ 
tical waveguides, (b) an array of circularly curved optical 
waveguides, (c) an array of circularly curved optical waveg¬ 
uides with a sinusoidal modulation, (d) Schematic diagram 
of the ultrafast laser inscription technique, (e) White-light 
transmission micrograph of the facet of an array of straight 
waveguides. Lattice constant a = 16 jjim. Each waveguide is 
single-mode at 780 nm wavelength. The central waveguide, 
indicated with the red circle, was excited for all measure¬ 
ments. 


trapped in optical lattices [33] , 

In this Letter, we use photonic lattices, fabricated us¬ 
ing the technique of ultrafast laser inscription , as a pow¬ 
erful platform to investigate the dynamics of W-S states 
in a strong static electric field and weak inter-site inter¬ 
action. We demonstrate strong localization of the W-S 
state, seen for the first time using curved photonic lat¬ 
tices, where the curvature is analogous to the inverse of 
a static electric field in the electronic case. When the 
electric field exceeds a threshold value, we observe that 
the W-S state becomes localised to a single lattice site. 
Importantly, we then also demonstrate that a strongly 
localised W-S state becomes delocalised when an appro¬ 
priate (specific frequencies and amplitudes) sinusoidal 
modulation is applied to the lattice. The latter result 
constitutes the first photonic-crystal analogue of photon- 
assisted tunneling, based on fabricated sinusoidal mod¬ 
ulations, offering a promising method for the generation 
of gauge fields in photonic lattices. 

The photonic lattice. The propagation of the electric 
field envelope in the material is governed by 




- + V{x' - xo{z),y) 
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where ^ depends on x',y and z, A = 27rA is the free- 
space wavelength and ^ V{x',y) = 

Vo{x'—Xn, y—ym) describes the refractive index mod¬ 
ulation in the transverse cross section where Vo{x',y) is 
the refractive index profile of a single waveguide at po¬ 
sition Xn^yn- The function xo{z) determines the trans¬ 
verse shift of the whole lattice depending on the propaga¬ 
tion distance By making a change of reference frame 


(a) 

. . Q 

(b) 

• 

(c) 

• • • 

(d) 

• 

(e) 

• 

(f) 

• 


FIG. 2: Intensity distribution at the output of 30-mm-long 
circularly curved ID lattices with radii of curvature (a) R = 
oo (i.e. straight lattice), (b) i? = 1.5 m, (c) R = 1.2 m, (d) 
R = 0.5 m, (e) R = 0.3 m and (f) R = 0.2 m. In each case, 
light was launched into the central waveguide. 


X = x' — xo{z) equation can be rewritten as 

$ (4) 
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with F = —nod‘^xo{z) and where ^{x^y^z) is the conse¬ 
quently transformed state [23l |34] . The transport of light 
in a circularly curved one-dimensional photonic lattice, 
with sinusoidal modulation (FigB c)) is then governed 
by the paraxial equation [23[ |34] 


8^ 

= 


—noAujQ sin(u;o^)^ 


( 5 ) 


where the lattice is bending along the x direction with a 
radius of curvature R. The amplitude and frequency of 
the ^-dependent ”ac” modulation are A and cjq respec¬ 
tively, no is the refractive index of the substrate mate¬ 
rial and An(x, y) is the transverse refractive index pro¬ 
file. Eq. (§ is analogous to the Schrodinger equation 
of a particle with effective mass no and charge e mov¬ 
ing in a ID periodic potential V{x) = — An(x), with an 
external (artificial) electric field E = Edc + where 
cE = uq/RF no^co’o sm{uJoz). Here, 2 ; plays the role of 
time, cEdc = and cEac = uqAuJq sm{ujQz) (see Eqs. 

Q and 0). 

For well-confined single-mode waveguides, Eq. can 
be solved using the tight-binding approximation. When 
only the lowest band is excited, Eq. gives the coupled¬ 
mode equations 

d^s 

i^r— = + ^>*- 1 ) -s{a + K sin(a;o2:))$s, (6) 

oz 

where is the electric field amplitude in the 5:th waveg¬ 
uide, n is the nearest-neighbour coupling constant along 
the x-axis, a = n^a/RX is the strength of the external 
linear potential and K = n^Auj^a/X is the strength of the 
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FIG. 3: (a) Micrograph of the facet of a square lattice. For all 
measurements, light was launched into the circled waveguide. 
(b)-(g) Intensity distribution at the output of 10-mm-long 
circularly curved square lattice with radii of curvature (b) 
R = oo (i.e. straight lattice), (c) i? = 1 m, (d) R = 0.7 m, 
(e) R = 0.5 m, (f) = 0.4 m and (g) R = 0.3 m. tvx — 0.085, 

mm“^, hiy = 0.095 mm“^, hiN = 0.019 mm“^, a=15 /xm. 

ac modulation. When only the dc field is present, and 
2K.RX/noa <C 1, the state excited at the input becomes 
localised to one site. When adding the ac field, this lo¬ 
calised state can interact with its neighbouring sites when 
the resonance condition 

uqu = a = rioa/RX^ z/G Z, (7) 

is satisfied. This leads to an effective tunneling given by 

ma 

( 8 ) 

K \UJ0 J 

where Jy is the Bessel function of order v. 

Wannier-Stark localisation To investigate Wannier- 
Stark localisation, fifteen ID lattices (lattice constant a = 
16 /xm) were fabricated using ultrafast laser inscription, 
Fig. [^d). In these lattices, no modulation was created 
(M = 0), however the radius of curvature of the lattice 
was varied between 1.5 and 0.1 m (i7 = 1.5,1.4,..., 0.1 
m), see Fig. Bb). An additional straight lattice was 
also fabricated {R = oc), see Fig. [^a). The white- 
light transmission micrograph of the facet of a lattice is 
shown in Fig. [^e). The refractive index profile of each 
waveguide was controlled using the “slit-beam” shaping 
method [35]. Each waveguide was inscribed by translat¬ 
ing the 30-mm-long glass sample (Corning Eagle^^^^), at 



FIG. 4: Inverse participation ratio (IPR) is a measure of lo¬ 
calisation and is defined as the average of the absolute value 
of the fourth power of the wave function. The IPRs along 
the X and y directions, (IPR)cc and (IPR)^/ respectively, have 
been plotted as a function of the inverse of curvature radius R 
which is a measure of the strength of dc electric field. There is 
no effect of electric field along the y direction. Gomplete local¬ 
isation (IPR=1) is observed in the x direction as the electric 
field exceeds a threshold value. 

a translation speed of 8 mm-s“^, once through the focus 
of a 500 kHz train of 1030 nm femtosecond laser pulses. 
The laser inscription parameters were optimized to pro¬ 
duce waveguides that were single-mode and well confined 
at 780 nm. For a more detailed description of the waveg¬ 
uide fabrication procedure, see Ref. [36] , 

The nearest-neighbour coupling n was measured to be 
0.072 mm“^. Fig. shows the output intensity distri¬ 
bution measured for the lattices with radii of curvature 
R = oo (i.e. the straight lattice), 1.5 m, 1.2 m, 0.5 m, 
0.3 m and 0.2 m. It is clear from Fig. [^that the light be¬ 
comes increasingly localized as the radius of curvature is 
reduced, as would be expected from Eq. Q. To investi¬ 
gate this phenomenon further, we fabricated 10-mm-long 
2D lattices with a lattice constant a = 15 jam along both 
X and y axes, where each lattice curves only along the 
X direction. Eor these lattices, the measured coupling 
strengths along the x and y axes were 0.085 mm“^ {tix) 
and 0.095 mm“^ {ny) respectively. Erom simulations, the 
estimated value of next-nearest-neighbour coupling was 
0.019 mm“^ (/^at)- As can be seen from Eig|^ localisation 
occurs only along the x axis, the direction of the artificial 
electric field. To quantify localisation along the two axes, 
the inverse participation ratio (IPR) was calculated. The 
IPR is a measure of localisation and is defined as the in¬ 
verse of the absolute value of the average of the fourth 
power of the wave function. Eor our purpose, the IPR 
for the X axis was obtained by summing all the intensity 
values in each column to obtain a vector of values along 
the X axis. The IPR was then calculated using this vec¬ 
tor. The IPR along the y axis was calculated using the 
same procedure, but by summing the rows rather than 
columns. Eor a localised state, the IPR is equal to 1. As 
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FIG. 5: Intensity distribution at the output of circularly 
curved lattices with sinusoidal modulation. Spatial periods of 
oscillation are 3.9 mm (l:st column), 7.8 mm (2:nd column) 
and 11.7 mm (3:rd column). For each period the amplitude 
of oscillation, A, was varied as shown. See also Fig. §. 

can be seen from Fig. there is no effect of electric field 
along the y axis, as would be expected, but complete lo¬ 
calisation is observed along the x axis once the artificial 
electric field exceeds a threshold value. 

The localization phenomenon can also be explained us¬ 
ing the theory of waveguide optics. It can be shown [34] . 
using a conformal transformation, that a ID array of cir¬ 
cularly curved waveguides with periodic transverse re¬ 
fractive index profile is equivalent to an array of straight 
waveguides with a new refractive index profile. In the 
limit a/R 1, the new refractive index profile is the su¬ 
perposition of the original periodic index profile and a lin¬ 
ear ramp of refractive index, and the radius of curvature 
controls this ramp. In other words, the mode supported 
by each waveguide in the curved array has a different 
propagation constant p. As decreases, the difference 
in p increases, resulting in partial transfer of light into 
the nearest waveguides via evanescent coupling. After 
a threshold value of i7, there is no significant coupling 
between the nearest waveguides resulting in a complete 
spatial localization. 

Analog photon-assisted tunneling To observe the effect 
analog to photon-assisted tunneling, three sets of 30-mm- 
long lattices were fabricated with a sinusoidal modulation 
(Fig. 0: c)). For sets 1, 2 and 3, the periods zq = 27tIujq 
were set to 3.9 mm, 7.8 mm and 11.7 mm respectively, 
corresponding to = 1, 2 and 3 in Eq. 0 . For all sets, 
the radius of curvature and inter-site separation were set 
to R=120 mm and a=16 ym respectively. For each set, 
15 lattices were fabricated, and the amplitude of oscilla¬ 
tion, A, was varied; see Fig. The measured output 
intensity distributions are shown in Fig The ratio 
of the effective coupling /^eff and the coupling strength 
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FIG. 6 : Graphical representation of Fig. §. Variation of 
effective coupling with Kq = Kjojo. The solid lines are the 
absolute values of Bessel functions | Ji |, | J 2 1 and | J 3 l^respec- 
tively, which were predicted theoretically from Eq. 

for the array of straight waveguides k was determined 
by comparing the output intensity distribution for the 
modulated curved lattice (Fig. [^c)) and that for the 
straight lattice (Fig. [^a)). The effective tunneling is 
plotted graphically as a function of Kq = KfujQ in Fig[^ 
where it can be seen that the effective coupling rate has 
a characteristic (Bessel-function) dependency on Kq, as 
predicted by Eq. This is clear evidence that the tun¬ 
neling has been partially restored through an analogue 
of photon-assisted tunneling. As a final note, it should 
be stressed that significant tunneling was absent when zq 
was not an integer multiple of 3.9 mm. 

Conelusion In this Letter, we have demonstrated that 
an appropriately designed array of evanescently coupled 
curved optical waveguides can be used to observe a W-S 
state that is fully localized on a single lattice site. From 
the perspective of solid state physics, the localization is 
due to an analogue of a strong external dc electric field 
that breaks the degeneracy of the Bloch states and re¬ 
sults in a Wannier-Stark ladder. We also demonstrate 
that tunneling in such photonic lattices can be restored 
by applying an analogue of an ac electric field, and that 
the strength of this tunneling obeys a characteristic de¬ 
pendency on the frequency and amplitude of the ac mod¬ 
ulation, which is in excellent agreement with the exist¬ 
ing theory of photon-assisted tunneling. By further tun¬ 
ing the spatial dependence of the laser-fabricated mod¬ 
ulation, this method could be used to produce effective 
magnetic fluxes [MIS] in 2D photonic lattices. The in¬ 
terplay between such artificial fields and the presence of 
non-linearities opens a promising route for the study of 
interacting particles in large magnetic fields. 
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